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PLATES WITH INCOMPATIBLE PRESTRAIN OF HIGHER ORDER 


MARTA LEWICKA, ANNIE RAOULT, AND DIEGO RICCIOTTI 


Abstract. We study the effective elastic behaviour of the incompatibly prestrained thin 
plates, characterized by a Riemann metric G on the reference configuration. We assume 
that the prestrain is “weak”, i.e. it induces scaling of the incompatible elastic energy E h of 
order less than h 2 in terms of the plate’s thickness h. 

We essentially prove two results. First, we establish the E-limit of the scaled energies 
h~^E h and show that it consists of a von Karman-like energy, given in terms of the first order 
infinitesimal isometries and of the admissible strains on the surface isometrically immersing 
Gt 2 x 2 (i.e. the prestrain metric on the midplate) in R 3 . Second, we prove that in the scaling 
regime E h ~ with /3 > 2, there is no other limiting theory: if inf h~ 2 E h —» 0 then 
inf E h < C/i 4 , and if inf h~ 4 E h 0 then G is realizable and hence min E h = 0 for every h. 
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1. Introduction 

The purpose of this paper is to study effective elastic behaviour of the incompatibly pre¬ 
stressed thin plates VL h , characterized by a Riemann metric G given on their reference con¬ 
figuration. The incompatibility is measured through the energy E h given below (sometimes 
called the “non-Euclidean” elastic energy). 
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We will be concerned with the regime of curvatures of G which yields the incompatibility 
rate of order higher than h 2 , in plate’s thickness h. Indeed, in paper [6] we analyzed the 
scaling inf E h ~ h 2 and proved that it only occurs when the metric G 2 x 2 on the mid-plate 
can be isometrically immersed in M 3 with the regularity W 2,2 and when, at the same time, 
the three appropriate Riemann curvatures of G do not vanish identically (for details, see 
below). The relevant residual theory, obtained through T-convergence, yielded a bending 
Kirchhoff-like residual energy. 

In the present paper we assume that: 

(1.1) hr 2 inf E h -A 0 

and prove that the only nontrivial residual theory in this regime is a von Karman-like energy, 
valid when inf E h ~ h 4 . It further turns out that this scaling is automatically implied by (O) 
and inf E h / 0. Indeed, we show that if (11.111 then h 4 inf E h < C, and that h 4 inf E h —> 0 
if and only if G is immersible whereas trivially rriin E h = 0 for all h. 

This scale separation is contrary to the findings of [29] valid in the Euclidean case of 
G = M3, where the residual energies are driven by presence of applied forces f h ~ h a . In 
that context, three distinct limiting theories have been obtained for E h h,P with (3 > 2 
(equivalently a > 2). Namely: f3 € (2,4) corresponded to the linearized Kirchhoff (nonlinear 
bending) model subject to a nonlinear constraint on the displacements, (3 = 4 to the classical 
von Karrnan model, and (3 > 4 to the linear elasticity. The present results are also contrary 
to the higher order hierarchy of scalings and of the resulting elastic theories of shells, as 
derived by an asymptotic calculus in [45] . The difference is due to the fact that while the 
magnitude of external forces is adjustable at will, it seems not to be the case for the interior 
mechanism of a given metric G which does not depend on h. In fact, it is the curvature 
tensor of G that induces the nontrivial stresses in the thin film. The Riemann tensor of 
a three-dimensional metric has only six independent components, namely the six sectional 
curvatures created out of the three principal directions, which further fall into two categories: 
including or excluding the thin direction variable. The simultaneous vanishing of curvatures 
in each of such categories correspond to the two scenarios at hand in terms of the scaling of 
the residual energy. 

1.1. Some background in dimension reduction for thin structures. Early attempts 
for replacing the three-dimensional model of a thin elastic structure with planar mid-surface 
at rest, by a two-dimensional model, were based on a priori simplifying assumptions on the 
deformations and on the stresses. Later, the natural idea of using the thickness as a small 
parameter and of establishing a limit model was largely explored; we refer in particular to 
the works by Ciarlet and Destuynder who set the method in the appropriate framework of 
the weak formulation of the boundary value problems na sa, proved convergence to the 
linear plate model [23] in the context of linearized elasticity, and obtained formally the von 
Karrnan plate model from finite elasticity [9] . See also [55l [56 ] for the time-dependent case 
and convergence results, and [10] for a comprehensive list of references. 

The issue of deriving two-dimensional models valid for large deformations, by means of an 
asymptotic formalism, was subsequently tackled by Fox, Sirno and the second author in |26| . 
They showed, in the context of the Saint Venant-Kirchhoff materials subject to appropriate 
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boundary conditions, how to recover a hierarchy of four models. This hierarchy, driven by 
the order of magnitude of the applied loads, consisted of: the nonlinear membrane model, the 
inextensional bending model, the von Karman model and the linear plate model. The models 
thus obtained still required a justification through rigorous convergence results. In [34], Le 
Dret and the second author used the variational point of view and proved T-convergence 
of the three-dimensional elastic energies to a nonlinear membrane energy, valid for loads of 
magnitude of order 1. We remark that the expression of the limiting stored energy therein 
consisted of quasiconvexification of the three-dimensional energy, first minimized with respect 
to the normal stretches. This allowed to recover the degeneracy under compression; a feature 
that is otherwise missed by formal expansions. We further mention that for 3d—Id reduction 
a similar point of view had been introduced by Acerbi, Buttazzo and Percivale in [1]. 

A key-point for deriving rigorously the above mentioned nonlinear bending model has been 
the geometric rigidity result due to Friesecke, James and Muller [28]. In the similar spirit, the 
same authors justified the von Karman model, the linear model [29] and also they introduced 
novel intermediate models, in particular in the range of energies - or equivalently of loadings 
between the scaling responsible for bending (/3 = 2) and the von Karman scaling (/3 = 4). 
In this range of models, the three-dimensional stored energy appears in the limit stored 
energy through its second derivative at rest. Scaling the energy with exponents /3 other than 
integers had been explored for the membrane to bending range in m leading to convergence 
results for 0 < /3 < 5/3 while the regime 5/3 < /3 < 2 remains open and is conjectured to 
be relevant for crumpling of elastic sheets. Other significant extensions concern derivation of 
limit theories for incompressible materials pH m E3 S3, for heterogeneous materials m , 
through establishing convergence of equilibria rather than strict minimizers [52H62L [3611531 137] . 
and finally for shallow shells [39 ]. 

Extension of the above variational method valid in the framework of the large deformation 
model was conducted in parallel for slender structures whose midsurface at rest is non-planar. 
The first result by the second author and Le Dret [35] relates to scaling (3 = 0 and models 
membrane shells: the limit stored energy depends then only on the stretching and shearing 
produced by the deformation on the midsurface. Another study is due to Friesecke, James, 
Mora and Muller [27] who analyzed the case (3 = 2. This scaling corresponds to a flexural shell 
model, where the only admissible deformations are those preserving the midsurface metric. 
The limit energy depends then on the change of curvature produced by the deformation. 
Further, the first author, Mora and Pakzad derived the relevant linear theories (/3 > 4) 
and the von Karman-like theories (/3 = 4) in [41], and subsequently proceeded to finalize 
the analysis for elliptic shells in the full regime [3 > 2 in 02 ]. A similar analysis has been 
performed in case of the developable shells in [32] leading to the proof of the collapse of all 
residual theories to the linear theory when f3 > 2. Following these findings, a conjecture was 
made in 05] about the infinite hierarchy of shell models and the various possible limiting 
scenarios differentiatied by rigidity properties of shells. Let us recall that a comprehensive 
body of work had been previously devoted to the asymptotic derivation of shell models in 
the small displacement regime under clear hypotheses on the model taken for granted, three- 
dimensional or already two-dimensional and containing the thickness as a parameter. Several 
models were recovered by Ciarlet and coauthors 0H |15j [16], by Destuynder [221 M and by 
Sanchez-Palencia and coauthors [59] EM El [50]. Sanchez-Palencia, in particular, theorized 
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the role and interplay of the midsurface geometry and of the boundary conditions [58], as 
well as underlined the singular perturbation behavior. We refer to HH for many additional 
references. 

Most recently, there has been a sustained interest in studying similar problems where 
the shape formation is not driven by exterior forces but rather by the internal prestrain 
caused by e.g. growth, swelling, shrinkage or plasticity [331 EH [62]. Variants of a thin plate 
theory can be used to study the self-similar structures which form due to variations in an 
intrinsic metric of a surface that is asymptotically flat at infinity [ 2 ], and also in the case of 
a circular disk with edge-localized growth [25] , or in the shape of a long leaf [48] . Ben Arnar 
and coauthors formally derived a variant of the Foppl-von Karman equilibrium equations 
from finite incompressible elasticity [201 El]: they use the multiplicative decomposition of 
the gradient proposed in [57] similar to ours. They take cockling of paper, grass blades and 
sympatelous flowers as examples mm- 

A systematic study of the possible limit problems when a target metric is prescribed was 
undertaken by the first author and collaborators: a residually strained version of the Kirchhoff 
theory for plates was, for the first time, rigorously derived in [54] under the assumption that 
the target metric is independent of thickness. This analysis was completed in [ 6 ], resulting 
in a necessary and sufficient condition that the elastic prestrained energy scales as h?. The 
object of the present paper is to study higher order prestrains. 

Let us also mention that in [39j |40l 03] similar derivations were carried out under a different 
assumption on the asymptotic behavior of the prescribed metric, which also implied energy 
scaling fr in different regimes of /3 > 2. In [39] it was shown that the resulting equations 
are identical to those postulated to account for the effects of growth in elastic plates @ 8 ] and 
used to describe the shape of a long leaf. In [43] a model with a Monge-Ampere constraint 
was derived and analysed from various aspects. Other results concerning the energy scaling 
for the materials with residual strain are derived in [4], where by imposing suitable boundary 
data, conditions of mm are not satisfied and hence the residual energy scales larger than 
h? . depending on the type of these boundary data (see also [62]). 

1.2. The set-up and notation. Let fl be an open, bounded, smooth and simply connected 
subset of R 2 . For 0 < h <C 1 we consider thin films Q h around the mid-plate fl: 

(1.2) Q h = [x = ( x',X 3 ); x' € X 3 € (—h/2,h/2)}. 

Let G : VL h —>• R 3x3 be a given smooth Riemann metric on Q h , uniform through the thickness: 

G(x',x 3 ) = G(x') for every {x 1 ,x^) € Q h , 

and let A = \[G denote the unique positive definite symmetric square root of G. Consider 
the energy functional E h : W 1 , 2 (Q h , R 3 ) —> R + defined as: 

(1.3) E h {u h ) = \ [ WiA/ybA- 1 ) dx. 

h Jn h 

The nonlinear elastic energy density W : R 3x3 —>• R + is a Borel measurable function, assumed 
to be C 2 in a neighborhood of 50(3) and to satisfy, for every F € R 3x3 , every R € 50(3) 
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and with a uniform constant c > 0 , the conditions: 

(1.4) W(R) = 0, W{RF) = W(F), W{F) > cdist 2 (F, 50(3)). 

The first condition states that the energy of a rigid motion is 0, while the second is the 
frame invariance. They imply that DWil&z) = 0 and that 0 2 lT(Id 3 )( 5 , •) = 0 for all skew 
symmetric matrices 5 € so(3). The third assumption above reflects the quadratic growth 
of the density W away from the energy well 50(3). Note that these assumptions are not 
contradictory with the physical condition W(F) = oo for det F < 0. 

Throughout the paper, we use the following notation. Given a matrix F € R 3x3 , we 
denote its transpose by F f , its symmetric part by sym F = ^(F + F f ), and its skew part 
by skewF = F — sym F. By 50(n) = {i? € M” xn ; Ft = R~ l and det 5! = 1} we denote 
the group of special rotations, while so(n) = {T € M nxn ; sym F = 0} is the space of skew- 
symmetric matrices. We use the matrix norm |T| = (trac e(F t F)) 1 / 2 , which is induced by the 
inner product (T\ : F 2 ) = trace(i ? ^i ? 2 ). The 2x2 principal minor of a matrix F € R 3x3 is 
denoted by ^2x2• Conversely, for a given 1*2x2 £ R 2x2 , the 3x3 matrix with principal minor 
equal 1*2x2 and all other entries equal to 0, is denoted by ( 1 * 2 x 2 )*- All limits are taken as 
the thickness parameter h vanishes, i.e. when h —>• 0. Finally, by C we denote any universal 
constant, independent of h. 

1.3. Some previous directly related results. It has been proved in [43] that: 

inf E h (u h ) = 0 
u h eVF 1 ’ 2 (fi h , R 3 ) 

if and only if the Riemann curvature tensor of G vanishes identically in fl h , i.e.: Riern(G) = 0, 
and when (equivalently) the infimum above is achieved through a smooth isometric immersion 
u h of the metric G on Fl h . Further, in [ 6 ] it is proved that: 

lim ^77 inf E h = 0 
0 h z 

if and only if the following Riemann curvatures of G vanish identically: 

(1.5) -R 1212 = -R 1213 = -R 1223 = 0 in Q h . 

More generally, the limit behavior of the rescaled energies h~ 2 E h has been investigated in [6] 
and it has been proved that their T-limit is given by the functional: 

My) = ^J n Q 2 A x '’^yyvb) dx r , 

effectively defined on the set of all y € W 2 , 2 (Fl, M 3 ) such that (Vy)*Vy = G 2 X 2 - The quadratic 
forms Q. 2 ,a{x' 1 ') are given by means of the energy density W as in (13.41) . The Cosserat vector 
b € W 1,2 D 3 ) is uniquely determined from the isometric immersion y by: 

( 1 . 6 ) Q l Q = G where Qe 1 = d\y, Qe 2 = 82 y, Qe 3 = b, with detQ > 0. 

Observe that the functional I 2 is a Kirchhoff-like fully nonlinear bending energy, which in 
case of Ge 3 = e 3 reduces to the classical bending content quantifying the second fundamental 
form (Vy)*V 6 = (Vy)*VlV on the deformed surface y(f2) with the unit normal vector N. 
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We recall that by Theorems 5.3, 5.5 and Corollary 5.4 in [6], the negation of condition 
((1.511 is equivalent to minl 2 > 0. For this reason, 111.5)1 is equivalent to the existence of a, 
necessarily unique and smooth, isometric immersion yo : ^ —>• R 3 of G 2 X 2 , such that: 


(Vyo)*V?/o — G*2 x2 
sym((Vy 0 ) i V6 0 ) = 0 . 


Above, the smooth vector field bo and the smooth matrix field Q o are given as in (11.61) : 

(1.8) QoQo = G, Qoei = diy 0 , Q 0 e 2 = d 2 y 0 and Q 0 e 3 = b 0 with det Q 0 > 0. 
Equivalently, denoting the inverse matrix G _1 = [Gr ij ']jj:i.. 3 , we have: 

(1-9) b 0 = --^(G 13 d m + G 23 d 2 yo) + -±=N. 

Uniqueness of the immersion yo hi (UUP follows from Theorem 5.3 in |6| which shows that 
the second fundamental form of the surface yo(fl) is given in terms of G. Therefore, both 
fundamental forms are known. Also, the second equation in (11,71) comes from the fact that 
the kernel of each quadratic form Q 2) a consists of so( 2). 


1.4. New results of this work. In view of the above statements, in this paper we investi¬ 
gate smaller energy scalings and the limiting behaviour of the minimizing configurations to 
E h under condition (USD- We first prove (in Lemma 12.11) that (USD , which as we recall is 
equivalent to (11.11) . implies: 

inf E h < Ch A . 

We then derive (in Theorem IQ and Theorem 14.11) the T-limit I 4 of the rescaled energies 
h~ 4 E h , together with their compactness properties. 

Namely, let yo be the unique immersion satisfying (US), where bo is ns in (1 1 .8 [) . Let 
do : —>• R 3 be the smooth vector field given in terms of yo by: 

(1.10) {Qodo, ei) = — {d\b 0 , bo), (Qodo,e 2 ) =—(d 2 b 0 ,bo), (Qodo,e3)=0. 

The limit Z 4 is then the following energy functional: 

MV, S) = ^J Q 2 ,a (®',S+i(W)*W + ^(V6o)*V6o) dx' 

(1-11) + ^L Q2 ’ A { - Vy °^ V P + ( Vy ) tv ^o) dx' 

+ t^oL Q2,a ( x ' : (vy ° )<vtr ° +dx/ ’ 

acting on the space of finite strains: 

§ € cl L 2 {sym((Vyo) < Vt(;); w € hF 1 , 2 (U,M 3 )} 
and the space of first order infinitesimal isometries: 

V G W 2,2 (f2,R 3 ) such that: sym((Vyo)*VU) 2 x 2 = 0, 
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where the vector field p £ VF 1 , 2 (fl,R 3 ) is uniquely associated with each V by: (Vyo)V = 
-(W )*& 0 and (b 0 , p) = 0 . 

The spaces consisting of § and V contain the information about the admissible error 
displacements, relative to the leading order immersion yo, under the energy scaling E h ~ fi 4 . 
We discuss the geometrical significance of V and § and of various bending and stretching 
tensors in the first two terms of T±(V, §) in section [5j We further prove in Theorem 16.21 that 
the last term in HI.Ill) , which is obviously constant and as such does not play a role in the 
minimization process, is precisely given by the only potentially nonzero (in view of (11.5)) ) 
curvatures of G , namely: 

sym ((Vyo)*VcT 0 ) + (V6o)*V6 0 = 

We may thus write, informally: 

Z/i{V, S) = - / Q 2 ,a{x', stretching of order h 2 ) da/ + — / Q, 2 ,a(x', bending of order h) da/ 

2 Jn 24 Jn 

H-—— / Q 2 a(x', Riemann curvature of G) da/. 

1440 J n ’ 

In particular, since all three terms above are nonnegative, this directly implies that the 
condition lim/^o inf E b =0, which is equivalent to minX 4 = 0, is further equivalent to the 
immersability G, i.e. the vanishing of all its Riemann curvatures Riem(G) = 0 in Q h . 

1.5. Acknowledgments. M.L. was partially supported by the NSF grant DMS-0846996 
and the NSF grant DMS-1406730. A part of this work has been carried out while the first 
author visited the second author at the Universite Paris Descartes, whose support and warm 
hospitality are gratefully acknowledged. 
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2. The scaling and approximation lemmas 
We first introduce the following notation. Let Bq(x') be the matrix field satisfying: 

( 2 . 1 ) B^ei = d\bo, B$e2 = an<4 -^ 0^3 = do, 

where do is given by (|1.10l) . Observe that in this way QqBq is skew symmetric. Indeed, it 
has the following block form: 


( 2 . 2 ) 


QlB 0 


’ (Vyo)'V6 0 

(Vyo)* do 

(bofVbo 

(bo, do) 


and by (11.71) we see that (Vyo)* V 60 £ so( 2 ) is skew symmetric, while by (jl.lOD : (VyoY do = 
-(S/boYbo and ( 6 0 ,d 0 ) = 0 . 

Lemma 2.1. Condition 11.51) implies: inf E h < Ch A . 

w^ 2 (p. h , R 3 ) 


Proof. Let us construct a sequence u h € W 1 , 2 (£l h , R 3 ) that has low energy. Let: 

2 

(2.3) u h (x', x 3 ) = y 0 (x') + x 3 b 0 {x') + ydo (x 1 ), 
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in fact each u h is the restriction on its domain of the same deformation. We have: 

x 2 

Wu h (x',x 3 ) = Qo(x') + x 3 B 0 (x') + -^Dq(x'), 
where the matrix field Dq(x') € M 3x3 is given through: 

D 0 (x')ei = did 0 , D 0 (x')e 2 = d 2 d 0 , D 0 (x')e 3 = 0, 

so that: 

Vu h A~ l = Q 0 A- l + x 3 BqA ~ 1 + ^DqAT 1 . 

For brevity, denote F h = S7u h A~ 1 . Obviously, F h decomposes as: 

(2.4) F h (x', x 3 ) = QoA~ 1 {x'){J.d 3 + x 3 5(x') + x§T(x')) = (Q 0 A^ 1 (x , ))G h (x', x 3 ) 

with S = A^QqBqA " 1 , T = i A~ 1 QqDoA ~ 1 and G h = M 3 + x 3 S + x\T. Since QoA ~ 1 € 
50(3) by construction, frame indifference implies that W(F h ) = W (((G h ) t G h ) 1 ^ 2 ). Note 
that 5 is skew symmetric, because QqBo is skew symmetric. Therefore, ( G h YG h and the 
expansion of its square root do not contain terms linear in X 3 . Indeed, letting I\ = T+T t — S 2 \ 

((G h ) t G h )(x l , x 3 ) = Id 3 + 4K{x') + 0{x |) 

and: 

((OVo^/V, x 3 ) = Ids + y K ^ + 0 ( x 3)- 

As a consequence, using W (M 3 ) = 0 and DW (M 3 ) = 0, we obtain: 

4 

W(F h ) = W ^((G h ) t G h ) 1 ^ = ^D 2 W(Id 3 )(K,K) + 0(x 5 3 ). 

Using m, we get 

E h (u h ) = \ f W(F h ) dx < Ch 4 , 
h Jn h 

which accomplishes the proof of the lemma. ■ 


In Lemma 12.11 above, we constructed deformations whose gradient was sufficiently close 
to Q 0 + x 3 Bq, to provide the energy of the order h 4 . Conversely, in Corollary 12.31 below, 
we establish that the gradients of deformations u h whose energy scales like / 1 4 , are close 
to Q 0 + x 3 Bq modulo local multiplications by R h \x') € 50(3). Corollary 12.31 makes this 
statement precise and gives an estimation on \7R h as well. 

For any V which is an open subset of 0, we let V h = V x (— h/2,h/2) and we define the 
local energy functional by: 

E h (u h ,V h ) = l [ W(yu h A -1 ) dx. 
h J V h 
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Lemma 2.2. Assume There exists a constant C > 0 with the following property. For 

any u h € W 1 , 2 (V h , M 3 ), there exists R h € 50(3) such that: 


(2.5) ^ J \s7u h (x ) - R h (Q 0 {x) + x 3 B 0 (x')) f dx < C (V 1 ^, V h ) + . 


T/ie constant C is uniform for all V h which are bi-Lipschitz equivalent with controlled Lips- 
chitz constants. 


Proof. By assumption (11.41) . we have: 

(2.6) E h (u h ,V h ) >^J dist 2 (WOT 1 ,50(3)) dx. 

This suggests performing a change of variables in order to use the nonlinear geometric rigidity 
estimate [29]. For any u h € W 1 , 2 (V h , M 3 ), we let v h = u h o Y _1 with Y : V h —>• Y(V h ) = 
U h C M 3 given as in (12.31) . namely: 

2 

T (x', x 3 ) = yo(x') + x 3 b 0 (x') + y cfo(x'). 

Obviously, v h G W 1,2 (W ,l ,M 3 ) and: 

(2.7) V/r^x^xs) =V/(z)(VhA- 1 )( a ;',i3), *:=Y(x',x 3 ). 

Let 5' = 5 0 Qo 1 and T' = \D 0 Q~ l . Note that S' = B 0 Q~ l = (Q^Q- 1 ) = -Q" 1 -*#* 

in view of QqBq £ so(3). Therefore S' € so(3). Computations as in Lemma 12.II now give: 

2 

(2.8) VY(x', x 3 ) = Qo(x') + x 3 B 0 {x') + ^O 0 (x'), 
and: 

VYA” 1 = (Id 3 + x 3 5'(x') + x|T'(x')) (Qo^T 1 ). 

We see that for h small, det(VYA _1 ) > 0. Further, the left polar decomposition VYA -1 = 
(VYAl _1 (VY^4 _1 ) t ) 1 ^“ R, allows us to write: 

VYA” 1 = (Id 3 + x 2 M(x', x 3 ))R(x', x 3 ), 

where M = 0(1) is a symmetric matrix held and R £ SO( 3). Again, the symmetric term 
does not contain any term linear in x 3 . Therefore: 

dist (Vv h VYA-\ 50(3)) = dist (Vv h (ld 3 + x| M)R, 50(3)) 

= dist ^Vu fc (Id 3 + x\M), 50(3)) > cdist 50(3)(Id 3 + x§M)" 1 ) 

> cdist (Vv h , 50(3)) +0(xl). 

Now, let J = jdetVY o Y _1 | 1 . By (12.7(1 and the above computation: 

[ dist 2 (\/u h A~ 1 , 50(3)) dx > c [ dist 2 ('Vv h , 50(3)) J dz — c [ x\ dx. 

Jv h ' ' Ju h ' ' Jv h 
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In other words, since J > c > 0: 


1 

h 



50(3)) dx + h 3 \V h \ > 



dist 2 (Vv h , 50(3)) 


dz. 


By [29] . there exists C > 0 with the following property. For any v h € W 1,2 (U h , M 3 ), there 
exists R h € 50(3) such that: 


O [ dist 2 (\7v h , 50(3)) dz > f 
Ju h v ' Ju h 


Vv h - R h 


2 

dz. 


The constant C can be chosen uniformly for domains U h which are bi-Lipschitz equivalent 
with controlled Lipschitz constants. By (12.61) and the reverse change of variables which 
satisfies J -1 > c > 0 and |VY| < O, we obtain: 


C ( E h (u h ,V h ) + h 3 \V h \) > f 


1 

h J\?h 


\7u h - R h X7Y 


dx, 


again with a constant C uniform for domains V h that are bi-Lipschitz equivalent with con¬ 
trolled Lipschitz constants. This accomplishes the proof in view of (|2.8D . ■ 


Corollary 2.3. Assume 11.51) and let u h be a sequence of deformations such that: 

lim h~ 2 E h (u h ) = 0. 
h ->-0 

Then, there exist matrix fields R h € W 1,2 (VL, 50(3)) such that: 

1 


h 


J Vu h (x) - R h (x') (Q 0 (x') + x 3 B 0 (x')) 2 dx<C ( E h (u h ) + hf 


(2.9) 

and: 

( 2 . 10 ) 


The proof follows the lines of [ '291 [TL 55] . with necessary modifications in view of the 
expected error of the order /i 4 . For completeness, we will present the details in the Appendix. 


VR h (x') ' dx' < ^ (E h (u h ) + /i 4 ) 


3. The lower bound 

Theorem 3.1. Let u h € W 1,2 (Ll h , M 3 ) be a sequence of deformations satisfying E h {u h ) < 
O/i 4 . Then there exists a sequence of translations c h €E M 3 and rotations R h € 50(3) such 
that the associated renormalizations: 

(3.1) y h (x',x 3 ) = (R h ) t u h (x',hx 3 )-c h € 1T 1,2 (9 1 ,M 3 ) 

have the following properties, where j/o and bo are the unique solution to (fT?p (fOp. All 
convergences hold up to a subsequence: 

(i) y h —»■ yo in lF 1 , 2 (fl 1 ,M 3 ) and \d 3 y h —>■ bo in L 2 (r 2 4 ,R 3 ); 
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(ii) the scaled average displacements: 

_ i 

(3.2) V h (x') = ri (y h (x',x 3 ) - ( y 0 {x ') + hx 3 b 0 (x'))^ dx 3 

" 2 

converge in W 1,2 (Q, M 3 ) to a limiting field V G W 2,2 (Q, M 3 ), satisfying the constraint: 

(3.3) sym ((Vy 0 )*W) = 0; 

(iii) the scaled tangential strains: 

UymdVyoYVVb) 


h 


d2x2 


converge weakly in L 2 (fl,IR 2x2 ) to some § G L 2 {Pl,^ S y m) . 

(iv) Further, defining the quadratic forms Q 3 (F) = D 2 W (ld 3 )(F, F) and: 

(3.4) Q 2 ,a(x',F 2x2 ) = min{g 3 (^(zO _1 Mz , r 1 ); F G M 3x3 with F 2x2 = F 2x2 } , 

we have: 


lim inf ■ .. 
h->0 /i 4 


±E h (u h )>l 4 (V,S) = \f n Q2,A (x , ,§ + ±(VV) t VV + ^(Vb 0 ) t Vb 0 ' S J dx' 
(3-5) + 1 LJ n Q2 ’ A (Vyo ^P + ( W )' V ^>) dx' 

0,2,A (s', (Vyo)*Vd 0 + (Vf>o)*V6 0 ) dx', 


+ 


1440 


(3.6) 


where the vector field p G W 1,2 (g,M 3 ) is uniquely associated with V by: 

I (Vyo^P = -(W )‘& 0 
\ (bo , P) = 0 . 


Proof. 1. Corollary 12.31 yields existence of R h G W 1,2 (fl, ,SO(3)) such that (12.911 and (12.101) 
hold with Ch 4 and Ch 2 in their right hand sides, respectively. We rewrite these inequalities 
for the reader’s convenience: 


(3.7) 
and: 

(3.8) 


\ 


'n h 


V« h (x) - R h (x') (Qo(x') + x 3 B Q {x')) 


dx < Ch 4 


VR h (x') 


dx' < Ch 2 . 


To prove the claimed convergence properties for (13.11) . it is natural in view of (13.71) to set: 

R h = ^SO( 3 ) / Vu h (x)Q 0 (x , )~ 1 dx. 

Jn h 
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This projection is well defined, because for every x' € 0, in view of (13.71) : 


dist 2 (/ Vu ft Qo 1 dx, 50(3)) < / Vu h Q^ 1 dx 

VJo' 1 J Jn h 


- R h {x') 


< C 

< C 


/ (yu h Q^-R h ) dx +C / R h dx — R h (x') 

Jn h Jn h 


o' 1 


Vu h - R h (Q 0 + x 3j B 0 )) Qq 1 


dx + C 


i ? 1 


V) - / /?/ 

Jo 


<ci | W 1 - i? ft (Q 0 + x 3 B 0 )\ 2 dx + Cji?V) - / 
./o ' 1 Jo 


- T R 
o 


/i|2 


< C/i 4 + C\R n (x') - f R 

Jo 


h |2 


Now, taking the average on fl, by the Poincare-Wirtinger inequality and (13.81) . we get: 

dist 2 ( / Vu h Qo 1 dx, 50(3)) < O/i 4 + O [ \VR h \ 2 < Ch 2 , 

VJn' 1 / Jo 

which proves that the average f nh Vu^Qq 1 dx is close to 50(3) and that: 


(3.9) | / Vu h Qc 1 dx - i2 fc | 2 < Ch 2 . 

Jn h 

Moreover: 

-f \R h - R h \ 2 dx = \R h — R h \ 2 dx 
Jo J o ' 1 

(3.10) <C-f (\R h — -f R h \ 2 + \(-f R h ) — -/ V^Qo 1 ! 2 ) + / |J? h - / Vi/Q^l 2 

Jo ' 1 V Jo Jo Jo ' 1 / Jo h Jo ' 1 

<o/ |Vi? h | 2 dx + o/ |Vii h — R h (Qo + x 3 i?o)| 2 dx + Ch 2 < Ch 2 , 

Jo ' 1 Jo ' 1 

where the last estimate follows by (13.71) . (13.81) and (13.91) . 

Let now c h € M 3 be such that f n V h = 0 where V h is defined as in (13.21) . Denote by VhU h 
the matrix whose columns are given by d\y h , 82 y h and d 3 y h /h. Obviously: 

( 3 . 11 ) V, * 3 ) = (R h Y Vn ft (x', hx 3 ). 


Observe that: 


[ \y h y h -Qo\ 2 dx<ci \Wu h — R h Qo \ 2 dx 
Jo 1 Jo' 1 


<CU \Vu n - R tl {Q 0 + x 3 B 0 )\ 2 dx+f \x 3 R n B 0 \ 2 dx + f \R h - R h \ 2 dx) < Ch 


J O ' 1 J O ' 1 J O ' 1 

by (|3.7D and (13.101) . Therefore, V hy h converges in L 2 (D 4 ) to Qq. Further, the sequence 
{ y h } is bounded in 1T 1,2 (D 1 ), by the choice of c h . Passing to a subsequence we get that y h 
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converges weakly in W 1,2 (yi l ) and in view of the strong convergence of S7y h we have: 

1 


y —>• yo in W i,2 (fi,M d ) and —<9 3 y -A- bo in L 2 (fr,M 3 ). 


2. Note that, for every x' € Q: 


S7V h (x') = 


1 


1/2 

- 1/2 


Vhy h (x) - Q 0 (x')dx 3 


3x2 


(3.12) 


1 

h 


1/2 


- 1/2 


V h y h - ( R h YR h (Q 0 + hx 3 B 0 ) dx 3 j + ^ [((R h ) t R h - Id 3 )Q 0 ) 


= /i +4- 


The first term above converges to 0. Indeed: 

2 cr 


(3.13) 


4 L\n) < h 2 J n i 
C 


\(R h yVu h (x', hx 3 ) - {R h ) t R h {Q 0 {x') + hx 3 B 0 )| 2 dx 


3x2 


h\t r>hi 


3x2 


< 


h 2 


-f \Vu h (x',x 3 ) - R h (Qo + X 3 B 0 )\ 2 dx < Ch 2 . 
Jn h 


Towards estimating the second term in (|3.12|> . denote: 


S h = -((R h YR h -Id 3 ). 
h 


By (13.101) and (13.81) . it follows that: 




R h - R h 


C 


/ 


VR h I 2 < C. 


<C and \\VSTmn) < h2 

Passing to a subsequence, we can assume that: 

(3.14) S h S weakly in W 1 ' 2 ^), 
which implies: 

(3.15) I%^(SQ 0 ) 3x2 in L 2 (Q,M 3x2 ). 

Consequently, by (|3.12|) : 

(3.16) VV h ^(SQ 0 ) 3x2 in L 2 (0,E 3x2 ). 

As before, we conclude that V h converges in VP 1,2 (H) and that its limit V belongs to 
TT 2,2 (11,M 3 ), since W = ( SQq ) 3X 2 G IP 1,2 (S2). We now prove (13.31) . By definition of S h : 

(3.17) sym S h = ~{S h ) t S h , 

so in view of the boundedness of {S h } in VP 1,2 : 

||sym^|| L2(n) < C7i||S h ||i 4(n) < Ch\\S h \\ 2 wl , 2m < Ch. 

Consequently, S is askew symmetric field. But (Vyo)*VP = ( QqSQo) 2 x 2 , hence (13.31) follows. 
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For future use, let us define p G W 1 , 2 (f2, M 3 ) by: 

(3.18) [VF | p\ = SQ 0 . 

Since Qq[VF | p] G so(3), it is easily checked that p is given solely in terms of F by: 

f(Vy 0 )V=-(VF) t 6 0 

\(bo, P) = 0 . 


(3.19) 


3. We now want to establish convergence in (iii). In view of (13.1211 we write: 
1 sym (QoVF ft ) 2 x 2 (x / ) = ^ sym (Q£/{*) 2 x 2 + ^ sym (q^Qq^ 


(3.20) 


h 


2x2 


= 4 + 4 - 


We first deal with the sequence J 2 . By (I3.14|> . S h —>■ S in L 4 (fi) and so (|3.17l) implies: 

1 


h 


sym S h ^-^S t S=^S 2 in L 2 (Q). 


Therefore: 

(3.21) 


A^-\(Q t oS t SQo) 2 x 2 = -\AV) t VV in L\n). 


We now prove that 4 converges. Recall that by (13.201) . ()3.12l) and (13.111) : 


(3.22) 


h\t 


4 = synr (Qo/i ) 2X 2 = sym I Q 0 (R ) 


1/2 


- 1/2 


Z h (x',x 3 )dx 3 


2x2 


where the rescaled strains Z h are defined by: 

(3.23) Z h (x', x 3 ) = ( 4u h (x ', hx 3 ) - R h (x')(Q 0 (x') + hx 3 B 0 (x '))) . 

By (13.71) . the sequence {Z h } is bounded in L 2 (11 1 ,M 3 ). Therefore, up to a subsequence: 

(3.24) Z h —*• Z weakly in L 2 (11 1 ,M' 3 ). 

It follows that: 


(3.25) 


4 


J\ := sym ( Q^R) 1 Z(x',x 3 ) da ?3 J weakly in L 2 (fl). 

V J — 1/2 J 2x2 


which yields (iii) by (13.201) and (13.211) . 

4. We now aim at giving the structure of the weak limit § of ^ sym (QoVF ft ) 2x2 m terms 
of the limiting fields F and Z. We have just seen that: 


(3.26) 


s = Ji-^vvyvv, 


where J\ is given by (|3.25l) . As a tool, consider the difference quotients f s,h : 
f s ’ h (x',x 3 ) = (y h (x',x 3 + s) - y h {x',x 3 ) - hs (b 0 + h(x 3 + 
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and let us study for any s the convergence of f s,h as h —>• 0. In fact, we will show that 
f s,h p , weakly in W l, 2 {£l l , R 3 ). Write: 


<9 32 / (y, x 3 + t) - h(b 0 + h(x 3 + t)do) d t, 


/*-v,* 3 ) = p£ 

and observe that: 

jp ( <9 3 y /l - h(b 0 + hx 3 do)^j = i ^( J R h ) i Vu /l (x' ) /tx' 3 ) - (Q 0 + hx 3 B 0 )^ e 3 

= ^(R h y (vu h (x', hx 3 ) - i? h (Qo + hx 3 Bo)^J e 3 + S h (Q 0 + hx 3 B 0 )e 3 
= h(R h yZ h (x',x 3 )e 3 + S h (Q 0 + hx 3 B 0 )e 3 . 

The first term in the right hand side above converges to 0 in L 2 (n 1 ) because { Z h } is bounded 
in L 2 (n 1 ,M 3 ), while the second term converges to SQoe 3 = Sbo in L 2 (0 1 ) by (13.141) . Note 
that SQoe 3 = p by (|3.18]h Therefore, f s,h -A p in L 2 (H 1 ). 

We now deal with the derivatives of the studied sequence. Firstly: 

1 / 1 


<93 f s ’ h (x', x 3 ) =- (^ (< 93 y h (x', x 3 + s) - h(bo + h{x 3 + 
- ^ { d 3 y h ( x ', x 3 ) - h(b 0 + 
converges to 0 in L 2 (n 1 ). For i = 1,2, the in-plane derivatives read as: 


$/*’V, *3) = ^ /l(x 3 + s )) 


- (R h ydiU h (x' , /ix 3 ) - /is (+ h ( x 3 + ^ ) ) <9, 


1 


= - [{R h yz h {x', x 3 + s) - (^) t Z / *(x , ,s 3 )J ei 

1 / (5 A ) t i2 A (Qo + h(x 3 + s)B 0 ) - ( R h yR h (Q 0 + hx 3 B 0 )) e t 


+ 


h 2 s 

1 ( g 

~ [B 0 ei + h(x 3 + - 


The last two terms above can be written as: S h Boei — (x 3 + |) dido, hence by (13.241) : 

-r' 

s 


(9 if s ' h (x r , x 3 ) ^(Ry (z{x', x 3 T s) Z[x', x 3 )j ei 


+ SBoei — ( x 3 + -j dido weakly in 


where R € SO (3) is an accumulation point of the rotations R h . 

Consequently, f s,h —*• p weakly in W’ 1 , 2 (n 1 ,R 3 ) and, for i = 1 , 2 : 

(3.27) sdip = (. Ry(z{x',x 3 + s) - Z(x',x 3 )^ei + sSB 0 ei - s (x 3 + 0 dido, 
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which proves that Z(x', -)ej has polynomial form and that: 

(3.28) {&Z(x!,x, 3 )) 3x2 = (i? t Z(x',0)) 3x2 + x 3 (Vp- (SB 0 ) 3x2 ) + |W 0 . 

By (|3.24p . it follows that: 

Ji = sym 0)) 2x2 + sym (QoVd 0 ) 2 x 2 - 

With (|3.26l) . we hnally arrive at the following identity that links § and V and Z: 

(3.29) S(s') = sym (Q^Kf Z (a/, 0)) 2x2 + ^ sym (Q* Vd 0 ) 2x2 - 

5. We now prove the lower bound in (iv). Recall that by (13.231) : 

Vu h (x', hx 3 ) = R h (x')(Q 0 (x') + hx 3 B 0 (x')) + h 2 Z h (x',x 3 ). 

Since QqA ~ 1 € 50(3) we have: 

W{Vu h A~ l ) = W {[QQA~ 1 ) t (R h ) t Vu h A~~ 1 ^ = W(ld 3 + hj + h 2 Q h ), 

where: 

J(x',x 3 ) = x 3 A^ 1 {Q t 0 Bo)A~ l {x') € so(3), Q h [x\x 3 ) = A~ 1 Q t 0 (R h ) t Z h (x',x 3 )A~ 1 . 
Note that by (13.241) : 

G h (x',x 3 ) ~^G = A~ 1 Q t Q (R}‘)Z(x',x 3 )A~ l weakly in L 2 (n 1 ,M 3x3 ). 

Define the “good sets”: 

Glh = {a: € fi 1 ; h\G h \ < 1}. 

By the above, the characteristic functions converge to 1 in L 1 (D). Further, by frame 
invariance and Taylor expanding W on 

W(ld 3 + hj + h 2 g h ) = W{e- hJ { Ida + hj + h 2 Q h )) 

= W(Id 3 + h 2 {G h - X -J 2 ) + o(h 2 )) 

= \q 3 (h 2 {g h - l -j 2 )^+o{h*). 

Therefore: 

liminf yjE h (u h ) > liminf [ tn h W (ld 3 + hj + h 2 g h ) dec 

fe-K) h A h^o h A Jqi ' 

(3.30) = liminf i Qs(in, sym ( 5 ‘ - \j ‘)) dx 

-kl n M sym ( g -\ J i) Ax ’ 

by the weak sequential lower semi-continuity of the quadratic form Q 3 in L 2 and in view of: 
Ij-jh sym [g h — 2 ^ 2 ) ~symt? — -J 2 weakly in L 2 ^ 1 ). 
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Note that by (|3.18l) we have: (QqS , S 0 ) 2x2 = —(VV^V&o and that: 

J 2 = -J l J = -xlA-'BtBoA- 1 . 

Therefore, using (I3.28H . the right hand side of (|3.30|) is bounded below by: 

- J Q 2 A (x , , S ym^Q t o (R) t Z(x , , 0 ) + x 3 (Q t o Vp + (VV) t Vb 0 )+^(Q t 0 Vd o + (Vbo) t Vb o )) 2x2 ) dx 
= \j\ 22 ,a (x', I (a/) + x 3 III{x') + x§JI(a/)) dx. 

Above we used (13.291) and we denoted: 


(3.31) 


I(x') = 8 - ± sym ((Vyo^Vdo) + ^(VV^VV 

H(x f ) = ^sym ((Vy o yVdo) + ^(V&o^V&o 
III(x') = sym ((VyoYVp) + sym ((W^V&o)- 


Let C- 2 ,a(x') be the symmetric bilinear form generating the quadratic form Q 2 ,a{x'). Since 
the odd powers of x 3 integrate to 0 on the symmetric interval (— 1 / 2 , 1 / 2 ), we get: 

J 0.2,a(x', I(x 1 ) + x 3 III(x') + x§17(a/)) dx 

= f Q2,a(x' ,I(x')) dx'+ ( [ x\dx 3 ) f Q 2 ,a(x',III(x')) dx' 

Jn J-i /2 Jn 


A/2 


A/2 


+ (/ xj dx 3 ) Q2,a{x',II(x')) dx' + 2 ( / x\ dx 3 ) / C 2 ,a{x' , I(x'), II(x')) dx’ 


— 1/2 


- 1/2 


— / Q2,a(x', I) + — / Q2,a(x', III) + — / Q2,a(x', II) + — / C 2 ) a{x' , I, II) dx' 


12 


80 


12 


= j Q2 ,a(x',I+ -^II) dx' + j Q2 ,a{x',III ) dx'+-^ j Q 2 ,a(x',II) dx' 
= £i(V,S), 

by a direct calculation. This completes the proof of Theorem 13.II in view of ()3.30D . 


4. The upper bound 

We now complete the proof of I 4 being the T-limit of h^ i E h , by proving that the lower 
bound (|3.5j) is optimal. 

Theorem 4.1. Let V G !T 2 ’ 2 (fi,M 3 ) and § G L 2 (S7,M 2 */() satisfy: 
sym ((Vy 0 ) f V V) = 0, 

§GS :=cl L 2 {sym((Vy 0 )*Vu;); w G W ll 2 (fi,R 3 )}. 


(4.1) 
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Then there exists a sequence u h G M 3 ) such that assertions (i), (ii) and (in) of 

Theorem I.Ill are satisfied with R h = Id and c h = 0, and: 

(4.2) - 1 


limsup - jE h (u h ) < X 4 (V,S). 
h-> o h 


Proof. In the construction below, we will use the following notation. In view of (|3.4I) . for 
every F 2 X 2 € R 2x2 one can write: 


(4.3) 


Q 2 ,a(x', F 2x2 ) = mm |q 3 (A 1 {F 2x2 + sym(c®e 3 ))A 


where Ff x2 denotes the R 3x3 matrix whose principal 2x2 minor equals F 2x2 . We will denote 
by c(x’,F 2x2 ) the unique minimizer in (14.31) . Note that c(x',-) is a linear function of F 2x2 
and it depends only on its symmetric part (symi ? 2 x 2 )- 

1. Since § € S, there exists a sequence w h € W 1,2 (J1,R 3 ) such that: 


(4.4) 


sym ((Vy 0 )*V (w h + —do)) ->• § in L 2 (D, R 2x2 ) 


and without loss of generality we can assume that each w h is smooth up to the boundary of 
D, together with: 


(4.5) 


lim Vh\\w h \\ W 2 ,oo = 0. 
o 


Fix a small eo € (0,1) and let v h G VF 2,00 (D, R 3 ) be a sequence of Lipschitz deformations 
with the properties: 

v h ^V in VF 2,2 (fl,R 3 ), 

^'ll^ll VK 2 '°° — e 0 ; 

lim y-yl-JV G S4; v h (x') F(a/)}| = 0. 


(4.6) 


We refer to 

ph G w 1 ’ 00 ^, 

(4.7) 


and 

by: 


h—>-0 /l 2 

for the construction of such truncated sequence v h . Define now 


/ = (® 


t\-1 


-(Vv h yb 0 

0 


and also define the fields (f l G IF 1,00 (D,R 3 ), %q smooth and f h G L°°(D,R 3 ) such that: 
Qlt = ^c(x',2(Vy 0 ) t Vtv h + (Wv h YVv h ) - 


1 

T-o 

•+o 

> 

1_ 


■ (y v hyph - 

0 


i|^| 2 


Q^ko = c(x', (VyofVdo + (V&o)*V& 0 ) - 
Qof h = c(x', (VyoyVp' 1 + ( Vv h yVb 0 ) - 


( vb 0 yd 0 
\do \ 2 

(Vv h yd 0 
(fl 1 , do) 
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lim ||r* - r h \\ L 2 = 0, lim /h ll^llipi.- = 0. 

0 / i —>-0 


Finally, let r* € W 1 , 0 O (fi,K 3 ) be such that: 

(4.8) 

It follows from the definition of the minimizing map c, that: 

Q 3 (A~ 1 (2Q t 0 [Vw h | /] + [W* | /]*[v/ | /j^ 1 ) 

= Q 2) a (s', 2 (Vy 0 ) t Vw h + (V/)*V/) , 
Qs^A-^QHVdo | fc 0 ] + [V6 0 I dofiVbo | do])^' 1 ) 

= Q 2 ,a(V, (Vyo)*Vdo + (V6o)*V6o), 
Q 3 (4r 1 (2Q* [V/ | /] + 2 [V/ | /ftV&o | do] V" 1 ) 

= Q 2 ,a(V, (Vy 0 )*V/ + (V/)*V6 0 


(4.9) 


Moreover, we have the following pointwise bounds: 

|/1 < C|V/|, 

|V/| < C(\Wv h \ + |VV|), 

I/I < CdV^I + |Wf + |V/||/| + |/| 2 ) < C(|Vi/| + |V/| 2 ), 
|V/| < C(\Vw h \ + \V 2 w h \ + |V 2 /||V/| + |V/| 2 ). 


(4.10) 


2. Consider the sequence u h € bF 1,00 (f/M 3 ) defined as: 

2 

u h (x',x 3 ) = yo{x') + hv h (x') + h 2 w h (x') + x 3 b 0 (x') + ^-d 0 (x') 

+ ^fco/0 + fcc 3 /(x') + h 2 x 3 t\x') + ^rV). 
0 2 

For every (x',x 3 ) € fl 1 we write: 

V/(a/, hx 3 ) = Q 0 (x') + Zf (s', x 3 ) + Z$(x’, x 3 ), 

where: 


h 2 Xo, 


Z/x',x 3 ) = d[Vu h | /] + //Vu;' 1 | /] + hx 3 [Vb 0 | d 0 ] + —^-/Vd 0 | feo] + ^x 3 [V/ | 
Z 2 /x',x 3 ) = ^[v£ 0 I 0] + h 3 * 3 [V/ | 0] + ^[V/ | 0], 

o z 


Since Qo-4 1 E iSO^), we get: 


Vu ft A- 1 (s', hx 3 ) = QoA- 1 (ld 3 + A^Q^A- 1 + A^Q^A^ 
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and, in view of (14.61) . (14.81) and (I4.10p . there follows for h sufficiently small: 

|| A^QlZ^A - 1 + A~ 1 Q t 0 Z%A~ 1 \\ L c O 

< C (/i|| Vt>^||x,oo + h\\p ' |l°° + h 2 || Vu/* |l^oo + h 2 ||g*||Loo + h\\ V6 o||l°° + ^IMolli, 00 
+ h 2 ||Vdo||L°° + h 2 ||A:o||i,oo + h 2 \\ ’S/jf l ^L co + 11 11 z,°° + ^ 3 ||VA?o||i, 00 

+ /i 3 ||V^|| L oo + ^HVr^llioo) < Ce 0 . 

By the left polar decomposition, there exists a further rotation R € SO{ 3) such that: 

RS7u h A~ l = ((Id 3 + A- X Q\Z\A~ X + vl _1 QqZ^A^Y (Id 3 + A^Q^A' 1 + A" 1 A" 1 )) 

= (ld 3 + 2A" 1 symiQ^Z^A- 1 + A~ x (Z^f Z^ A~ x + 0{\Z £|)) ^ 

= Ids + sym (QiZ^A ” 1 + ^A~ x (Z^f Z* A~ x 

+ 0{\ sym (Q* Z*) + (Z'fZ'f) + 0(\Z%\). 


3. Consider the set: 

= {(x' ,x 3 ) € ft; v h (x') = V(x')}. 

Note that on ft/i we have: ff l = p and QqIVu^ 1 | j/'] G so(3). Using Taylor’s expansion, it 
follows that: 

A / W{Vu(x',hx 3 )A~ 1 )dx = ^- [ Q^A-^Q^ + ^zWz^A-^dx + Sl 
J \Lh J 

where the error term can be estimated by: 

|£f I < A f n |2sym {Q\Z\) + (Z^Zfl 3 + |Z 2 Y + |2sym (Q*Zf) + (Z^Z^Z^ dx. 
Now on U/j we also have, by (14.101) : 

|2sym(Q^Zf) + | < c(h 2 \ Vw h \ + h 2 \Wv h \ 2 + h 2 + h 2 \Vv h \ + h 2 |VV| + h 2 |r*|), 

|Z 2 ft | < Ch 3 (l + |Vg*| + |Vr*|) 

< C7i 3 (l + \Ww h \ + \V 2 w h \ + |vV*||Vv ft | + \Vv h \ 2 + |Vf*|), 
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and therefore, in view of (|4.5|h (14.81) . (14. tip and V £ W 2,2 : 


[ + d.x 

J £lh 

[ h 6 \Vw h \ 3 + h 6 \Vv h \ 6 + h 6 + h 6 \Vv h | 3 + h 6 |vV | 3 + h 6 |r *| 3 dx 
Jn h 


C_ 

F 

c 


< ^(Fiiv^iiicopFiiv^lUa^ + FllVFllia + Fini + FiiVKHis 


+ /i 6 ||vV|| L co ||VV ||| 2 + (Vhll^llLcc ) 3 


0 us h —y 0. 


Analogously: 

4 [ \Z%\ 2 dx<^j [ h 5 + (h\\Vv h \\ L °o) 2 h 4 \V 2 v h \ 2 + h 6 \Vv h \ i dx^0 as h -A 0, 
h Jn h h J nh 

i [ fisymiQDz' + iZffztWzZldx 

- F / (^ 5 |V^| 2 + F|vV*| 2 + F| W*| 2 + h 5 + /i 5 |VC| + h 5 |V 2 C| + /i 5 |f 
+ h 5 |Vi/| 2 |V 2 y| + h 5 |W||V 2 y| 2 ) dx < Ce 0 . 


We therefore conclude that: 


(4.11) 


limsup \£i \ < Ce q. 
h-*o 


4. Consider now the error due to integrating on the residual subdomain: 


£ 


h 

2 



dx < 


C_ 

F 


/ 


2 sym((5o^i) + (Zf)*Z{*| 2 + |Z 2 h | 2 dx. 


Observe that, since the matrix field [\7v h \ jf 1 ] is Lipschitz, we have: 


sym((5o[Vx /l | p*’'])(x / ) < C\\S7v h \\ w i,oa dist (x',{v h = V}) 

Ceo 


< 


h 


■dist ( x',{v h = C}) —>• 0 in L°°(Q). 


The last inequality above follows by a standard argument by contradiction. If there was a 
sequence x h € 12 such that dist(x /l , {v h = V}) > ch , this would imply that: |{x'; ir(x ') / 
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V(x')}\ > n B{x h , ch)\ > ch 2 , contradicting (|4.6j) . Consequently, by (I4.5jl . (|4.8jl . ([4.6D : 

1^2 I < tj / h 2 1 sym(Qo[Vu /l | p h ])| dx 

+ ^t [ h A \Vw h \ 2 + h A \Vv h \ A + h A \ VV| 2 + h A \f ^\ 2 + h A + h 6 \Vv h \ A dx 

h Jn 1 \n h 

< ^o(h 2 )\n A \ n h | + \^|i/2|| V ^|| i2 

+ c\n 1 \u h \\\Vv h \\ A L s + ch\\v 2 v h \\ L o 0 ^\\v 2 v h \\ L 2\n 1 \u h \ 1 / 2 + ^-{Vh\\f^\\ L o 0 ) 2 \n 1 \u h 
+ {h\\V 2 v h \\ L °°) 2 \\Vv h \\ 2 Li \n 1 \U h \ 1/2 ^0 as/i^0. 


Thus: 


lirnsup j^E h (u h ) < limsup i iQ 3 (.4 1 (sym(QoZ^) + A 4 ) dx + Ce 0 . 

h ^0 h >-0 J 

Now on fi/j we have: 

2sym(Q t 0 Z 1 /l ) + (ZffZ? 

2 

= 2/i 2 (sym(Qo[V w h | /]) + y sym(Qo[Vd 0 | fc 0 ]) + x 3 sym(Q(,[Vp | r*])) 

+ h 2 ([VC | pf [VC | p] + x 2 [Vb 0 | doHVfto | do] + 2 x 3 sym([VC | pf [V6 0 | do])) + S h , 


where the present error £ h is estimated by: 


|£ fc | < c(h 3 \VV\\Vw h \ + h 3 |VC| + /i 3 |VC||Vp1 + h 3 |VC||f*| 

+ h A \Vw h \ 2 + h 3 \Vw h \ + h A \Ww h \\Wp + /i 4 |Vu; h ||f*| + h 3 
+ h 3 |Vp| + h 3 |f*| + h A + /i 4 |Vp| + h A |f*| + /i 4 |Vpf + /i 4 |f*| 2 ) 
< C/» 2 (o(l)V^|VC| +eg|V 2 C| + o(l)Vh + o(l)e 0 Vhy 


(4.12) 
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Consequently: 

limsup —rE h (u h ) 
h->o h 

< limsup- [ Q 3 (sym(Qo[Vu> /i | g*]) + sym(Qo[Vd 0 | fa]) 

ft-j-o z JCl h v 1 

+ x 3 sym(Qo[Vp | r*]) + ^[VR | pf [VR | p] 

+ 7^3 [V&o I 4]*[V6o | cf 0 ] + x 3 sym([W | pf [V 6 0 | do]))-4 _1 ) dx + Ceo 

= limsup \ [ Q 3 (A- 1 (sym(Q t 0 [Vw h \ <?]) + hw | pf[W | p] 
h-+ 0 1 Jil h V 1 

+ 7^ x 3 sym(Qg[Vdo | fa]) + -x§[V 6 0 \ dof[V & 0 \ d 0 ])4 _1 ) 

+ Q 3 ( x 3 s ym(Qo[Vp | r' 1 ]) + x 3 sym([W | pf [V&o | do]))^" 1 ) dx + Ceo- 

Denoting by: 

h(x') = sym((Vy 0 f Vi/) + i(Vv ft )W, / 2 (®') = ^ sym((Vy 0 f Vd 0 ) + (V&ofV&o, 
we have: 

Q 3 (h _ 1 (/*(x') + sym(c(x',/i(x')) 0 > e 3 ) + xp^x') + x\ sym(c(x', J 2 (x')) 0 e 3 ))4 _1 ) 

= Q 3 (4 _ 1 ((/i(x') + xp 2 {x'))* + sym(c(x / , h(x') + xp 2 (x’)) <S> e 3 ))A^ 

= 0.2,A^ihix') + xp 2 (x'j^, 

where we have used the definition and linearity of the minimizing map c. Recalling the 
definitions of the curvature forms I(x'), II(x') and III(x') in (13.3111 . observe that I 2 {x') = 
2 II(x') and that p\ converges to I in L 2 by (|4.4jl . Hence: 

limsup j-^E h (u h ) < i J Q 2 ,a[i{x') + xpI(x'j S j dx + i J Q 2 : a{x 3 III(x ')) dx + Ceo 

= X 4 (D,§) + Ce 0 . 

Since cq > 0 was arbitrary, the proof is achieved by a diagonal argument. ■ 


5. Discussion of the von Karman-like functional (13.511 

Theorems no and 14.11 imply, as usual in the setting of T-convergence, convergence of 
almost-minimizers: 

Corollary 5.1. If u h € W 1 , 2 (0 h , M 3 ) is a minimizing sequence to h^’ i E h , that is: 

-L E \u h )-ini±E h ^=0, 


lim 

h->0 
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then the appropriate renormalizations y h = (R h yu h (x', hx%) — c h £ W 1 , 2 (fi 1 ,R 3 ) obey the con¬ 
vergence statements of Theorem 13.11 (%), (ii), (in). The convergence of h~ l sym((yyo) t \7V h ) 
to § in (in) is strong in L 2 (fl). Moreover, any limit (V, S) minimizes the functional Z 4 . 

Proof. The proof is standard. The only possibly nontrivial part is the strong convergence of 
the scaled tangential strains in (iii), which can be deduced as in Theorem 2.5 in m- ■ 

Let us now compare the functional (13.51) with the von-Karman theory of thin shells that 
has been derived in [IT]. Recall that when S is a smooth 2d surface in M 3 , the T-limit of the 
scaled elastic energies h _4 (^ f sh T / L(V'u ft ')) on thin shells S h with mid-surface S, is: 

(5.1) i 4 ,s(V,S) = l -j Q 2 (§ — \{A 2 ) tan )dy + ±-f Q 2 ((V(iJV) - AIl) tan )dy. 

J S J s 

Above, II stands for the shape operator of S and N is the unit normal vector to S. The 
subscript tan means taking the restriction of a quadratic form (or an operator) to the tangent 
space T y S. The arguments of 1\s are: 

(i) First order infinitesimal isometries V on S. These are vector fields V £ W 2,2 (S, M 3 ) 
with skew symmetric covariant derivative, so that one may define: 

(5.2) A £ W 1,2 (S, so(3)) with A(y)T = d r V(y ) \/y € S VT(=T y S; 

(ii) Finite strains § on S. These are tensor fields § € L 2 (S, M 2 ^) such that: 

(5.3) § = L 2 — lim sym(VWh)tan for some Wh € W 1,2 (S, M 3 ). 

h-t 0 

In the present setting, denote S = yo(£l) and observe that the 1-1 correspondence between 
V in (15. 2 p and V in (|3.3p is given by the change of variables V = Voy 0 . The skew-symmetric 
tensor field A on T y S is then uniquely given by: 

(5.4) A(y 0 (x'))d e y 0 = d e V(x') and Ab 0 =p Me € M 2 , 
and the finite strains in (15.31) are related to ([4.11) by: 

(S(y 0 (x'))d e y 0 , d e y 0 ) = (§(x')e, e) Me £ M 2 . 

Recall that the first of the two terms in the functional m measures the difference of 
order h 2 , between the (Euclidean) metric on S and the metric of the deformed surface. 
Indeed, the amount of stretching of S in the direction r £ T y S, induced by the deformation 
Uh = id + hV + h 2 w , has the expansion: 

\d T Uh \ 2 ~ \t \ 2 = h 2 {2{d T w,T) + \d T V\ 2 ^j +0(h 3 ) = 2h 2 |((symViD)r, t) - ^(A 2 r, r)^ +0(h 3 ). 

The leading order quantity in the right hand side above coincides with: 

((symVw)e, e) + ^(d e V,d e V) = (symVtc + ^(VR) t VR)e, e^, 

where we write r = d e yo, for any e £ M 2 . This is precisely the argument of the first term in 
24 (V,§), modulo the correction (V6o)*V6o (equal to the third fundamental form on S in case 
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&o = N), due to the incompatibility of the ambient Euclidean metric of S h with the given 
prestrain G on Q h . 

The second term in (15.11) measures the difference of order h, between the shape operator 
II on S and the shape operator I^' on the deformed surface (id + hV)(S) whose unit normal 
we denote by N h . The amount of bending of S, in the direction r € T y S , induced by the 
deformation Uh = id + hV can be estimated by m- 

(Id + /iA)- 1 n ft (Id + hA)r - nr = (Id + hA)- 1 (d T N h + 0(h 2 ))r - nr 

= (Id + hi )” 1 ((Id + hA)Ur + h(d r A)N + 0(/i 2 )) - nr 

= (Id - hA)h(d r A)N + Q(h 2 ) 

= h(d T A)N + Q(h 2 ) = h(y(AN) - in) + Q(h 2 ). 

The leading order term in this expansion coincides with the term (\7yo) t \7p+ (VE) f V 6 o when 
&o = N, because in view of (15.41) : 

(( d T A)b 0 ,T) = (( d e (Ab 0 ),d e y 0 ) - ((Ad e b 0 , d e y 0 ) = (( d e p,d e y 0 ) + ({d e b 0 , Ad e y 0 ) 

= ((Vyo^Vp e,e) - <(VE)*V6 0 e,e), 

where we again wrote r = d e yo € T yo ^ x i^S, for any e € M 2 . This is precisely the argument in 
the second term in Z 4 (V,§). 

In the next section we identify the geometric significance of the last term in (13.51) . 

6. The scaling optimality 
In this section, we prove the following crucial result: 

Theorem 6.1. Assume H1.5\) . together with: 

(6.1) sym ((Vyo^Vdo) + (VboYVbo = 0, 

where yo, bo and do are defined in \1. 7[ ). \1.8\) . \1.10\) . Then the metric G is flat, i.e. 
Riem(G ) = 0 in Q h . Equivalently: min E h = 0 for all h. 

Observe that when bo = N, then by ( 11 . 101 ) there must be do = 0, and hence condition (16.11) 
becomes: N = const. This is consistent with our previous observation that when Ge 3 = e 3 , 
then already condition (11.71) is enough to conclude immersability of G in R 3 . Equivalently, 
£* 2 x 2 is immersible in R 2 , so that indeed yo(Cl) must be planar in this case. 

Towards a proof of Theorem 16.11 recall that Riem(G) is the covariant Riemann curvature 
tensor, whose components 12 ... and their relation to the contravariant curvatures in R are: 

Riklm = 2 (dklGim + dimGkl — dkmGil ~ dilGkm ) + G np (T^T^ — T^ m T^) 

Riklm GisRfcim, 
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where we used the Einstein summation convention and the Christoffel symbols: 

( 6 . 2 ) VI = ±G ns (d k G sl + djG sk - d s G kl ). 

In view of the symmetries in Riem(G) of a 3-dimensional metric G , its flatness is equivalent 
to the vanishing of the following curvatures: 

-R1212, -Rl213) -Rl223, ^1313, -Rl323, ^2323- 
The proof of Theorem 16.11 is a consequence of the following observation. 

Theorem 6.2. Assume \1.5\ ) and let yo, bo and do be defined as in Ji. 7 | ), \1.10\) . Then: 

(6.3) sym (( VyoYVdo ) + (V6 0 )*V6 0 = 

Proof. 1. We have: 

R 1313 = —dnGss + G np (T n 13 T p 13 - T^I* ), 

i?2323 = ~\d22G 33 + G np (r% 3 r P 3 - T^T*,), 

Ri 323 = -\di 2G33 + G np (T n 13 T p 23 - Tf 2 T p 33 ). 

On the other hand, in view of (11.101) : 

Vi, j = 1,2 i ({diyo, djd 0 ) + (djyo, dido)') = ^ (dj(diy 0 , d 0 ) + dfidjyo, do)) - (<%2/o, do) 

= ~ifdijG33 — (dijyo, do) 

because: dj(dj;yo, do) + di(djyo, do) = — dij\bo\ 2 = —dijG 33 . Consequently, the formula (|6.3I) 
will follow, if we establish: 

(6.4) Vi, j = 1,2 {d ij yo,d 0 ) = G np T* j T* i and (dJo, d 3 b 0 ) = G np T? 3 T p j3 . 

2. Before proving (16.4p we gather some useful formulas. Note that diG = 2 sym ((diQfiQ) 
for i = 1,2. Therefore, by direct inspection: 

(6.5) Vi, j, k = 1,2 (<%y 0 ,4vo) = ^(diG kj + djG ki - d k G i3 ). 

Also, recall that condition (II.7|) is equivalent to (see [6j, proof of Theorem 5.3, formula (5.8)): 

(6.6) Vi, j = 1,2 (d i:j yo , b 0 ) = ^( diG j3 + djG i3 ). 

Therefore, for all i,j = 1,2: 

{djVo , dfio) = dfidjyo , b 0 ) - (<%y 0 , b 0 ) = \{diG j3 - djG i3 ), 

(dibo,b 0 ) = i diG 33 . 


Rl3l3 R\323 
Rl323 R2323 


(6.7) 
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We now express d^y o, dibo and do in the basis {8\yo, d 2 yo,bo}, writing: 

dijy 0 = ajjdiyo + a‘f j d 2 yo + a%b 0 , 

(6.8) dSo = P\ d\ yo + pfd 2 y 0 + pfb 0 , 

do = 7 1( 9iyo + 7 2 <922/0 + 7 3fo o- 
By (16.51) . (16.61) . (16.71) and (11.101) . it follows that: 

G { a ij^ a ij^ a ij) = G Qo ld ijyo = QlAfi/o 

= ^ + djGu — d\Gij, diG 2 j + djG 2 i — d 2 Gij,diGoj + djG^ , 

G(Pi,Pi,P?y = GQo'dibo = Qldbo = QAbo 

= - (diGi3 — diGi3,diG 2 3 — 8 2 Gi3,diG33 S j , 

G{ 7 1 ,7 2 ,7 3 ) t = GQo'do = Ql/o = (<9iG 33 ,3 2 G 33 , o)'. 

In view of (16.21) we then obtain, for all i, j = 1,2: 

(a^a^a%) = (r?,-, b?-, r?-), (PL$,Pf) = (r^r* 3 ,r^), ( 7 1 ,7 2 ,7 3 ) t = (r 33 >ri 3 ,r| 3 ), 

so that (16.81) becomes: 

dijyo = Ajd 12/0 + Aj 821/0 + r Ijbo, 

(6-9) 8ib 0 = ri 3 5 iyo + rld 2 y 0 + T%b 0 , 

do = B 33 3iyo + r| 3 a 2 y 0 + B 33 6 o- 

3. We now prove (16.41) . Keeping in mind that QqQo = G, the scalar products of expressions 
in (16.91) are: 


(dijy 0 , do) = (T^dnyo + r|6 0 , r p 33 d p y 0 + T 3 33 b 0 ) = G np T^T p 33 , 

(dibo , djbo) = {T^dnyo + Tlbo,V%d p yo + T%b 0 ) = G np A 3 T p j3 . 
exactly as claimed in (16.41) . This ends the proof of Theorem 16.21 and also of Theorem 16.1 


7. Two EXAMPLES 

In this section we compute the energy Z 4 (V,§) in the two particular cases of interest: 
G(x',X 3 ) = diag(l, 1, \(x')) and G(x' ,x 3 ) = A(x , )Id 3 . 

Let p be as in the definition (13.61) . Writing: p = oPdiyo + oPd 2 yo 4- a 3 bo, we obtain: 

G(a 1 ,a 2 ,a 3 ) t = -((8^, b 0 ),(d 2 V, bo), 0 )*. 

Consequently: 

(7.1) P = -G li (8iV, bo)8 lV o - G 2i (diV, bo)d 2 yo - G 3i (diV, b 0 )b 0 . 
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Lemma 7.1. Let A : Ll —>• R be smooth and strictly positive. Consider the metric of the form: 
G(x', X 3 ) = diag(l, 1, X(x')). Then: 

(i) G is immersible in R 3 if and only if: 

M\ = V 2 A — —— VA ® VA = 0 in LI, 

2 A 

while the condition M\ ^ 0 is equivalent to: ch 4 < inf E h < Ch 4 . 

(ii) The T -limit energy functional Z 4 in H3.5\) becomes: 


Vw € VL 1 ’ 2 (Q, R 2 ) Vu € W 2 ’ 2 (Ll, R) 

1 r 1 1 

T±{v,w) = - / Q 2 (symVw + -Vu <g> Vu + —-VA <g) VA) dx' 

2 ,/n 2 96A 

+ L/ n Q 2 (^) + 5 4 i / n S 2 (M A ) d x', 

where Q 2 is independent of x' and it is defined by < 2 2 ,m in 
Proof. Part (i) of the assertion has been shown in [ 6 j. For (ii), note first that: 


yo(x') = x' and Q 0 = A = diag(l, 1, VT). 


Consequently, directly from (13.41) we see that Q 2 ,a = 0 . 2 ,id , which we denote simply by Q 2 . 

Further, in view of (14.lj) . every admissible limiting strain § £ § has the form § = symVtc 
for some w € W 1 , 2 (f2,R 2 ). Also, without loss of generality, every admissible limiting dis¬ 
placement V is of the form: V = (0,0, v) for some v € lT 2 ' 2 (n,R). We now compute, using 
(USD, (E2|) and (1741) : 


b 0 = \/Ae 3 , 


Therefore: 


do — <9 2 A, 0), 


p = —y/\(div, d 2 v, 0 ). 


(VboYVbo 

(VyoYVp 


-^VA<g)VA, (VyoYVdo 

—^Vu <g> VA — V\V 2 v, 
2\/A 


= —-V 2 A, 
2 


(W) 4 V 6 0 


—-=Vu <g> VA. 

2\/A 


This ends the proof of Lemma 17.11 in view of (13.51) . 


Lemma 7.2. Let A : Ll — >• R be smooth and strictly positive. Consider the metric G(x',x^) = 
A(x / )Id 3 - Denote f = \ log A. Then: 

(i) Condition (EH) is equivalent to A/ = 0 , which is also equivalent to the immersability 
of the metric G 2x 2 in R 2 . 

(ii) Under condition m, condition (16.11) can be directly seen as equivalent to Ric(G) = 0 
and therefore to the immersability of G. 
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(iii) The T-limit energy functional in H3.5\) has the following form: 

W§)= ^ [ e- 2 / Q 2 (S+ hvv) t VV + ±-e 2 fVf®\7f) dx' 

^ Jn ^ 

+ 1 kJ n Q2 ( 2Ws ® V/ ~ V ~ (Ws ’ V/)Id2 ) dx ' 

+ 14ioX Q2(e/RiC(G)2x2 ' ) ^ 

where Q 2 is as in Lemma \l~l\ and where Ric(G ) 2x2 denotes the tangential part of the 
Ricci curvature tensor of G, i.e.: 


Ric(G) 2x2 — 


Rll R12 
Rl2 R22 


Proof. The part (i) has been deduced in [ 6 ], together with the expression: 

(7.2) Ric(G) = -(V 2 / - V/ < 8 ) V/)* - (A/ + |V/| 2 )Id 3 . 

We now consider the case when (jl.7l) holds. By (i) the metric G 2x 2 is immersible in M 2 and 
in particular N = e 3 . Writing V = (V\, V 2 , V 3 ), from (11.91) . (16.91) and (17.11) we obtain: 

b 0 = \/Ae 3 , d 0 = -(difdiyo + d 2 fd 2 y 0 ), P = --^=(diV 3 di y 0 + d 2 V 3 d 2 y 0 ). 

(V 6 0 ) 4 V 6 0 = e 2 / V/ <g> V/, (WyVbo = e f VV 3 ® V/. 

Further, observe that: dido = -{dufdiy 0 + d 2i fd 2 y 0 + difduyo + d 2 fd 2i y 0 ), and so: 

j{diyo , dido) = -j{Xd u f + ^diXdif + ^d 2 Xd 2 f) = ~(d u f + |V/| 2 ). 

In the same manner, we arrive at: 

-^(<9 2 y 0 1 dj 0 ) = ~{d 22 f + |V/| 2 ), j{d 2 y 0 , d 3 d 0 ) = - d 12 f, j(diy 0 , d 2 d 0 ) = -d 2l /. 

Consequently, (Vyo^Vefo is already a symmetric matrix, and: 

(VyoYVdo = —e 2 7(V 2 / + |V/| 2 Id 2 ). 

In particular, under condition A/ = 0, the formula (17.21) yields: 

sym(Vy 0 ) # Vdo + (V 6 0 ) f V 6 0 = e 2 / Ric(G) 2x2 , 


which we directly see to be equivalent with V/ = 0 and hence with Ric(G') = 0. This 
establishes ( ii ). 

We now compute the remaining quantities appearing in the expression of I 4 . Firstly: 

Vp = 7^^Vpo(VV3 ® VA) - -i=Vy 0 V 2 R 3 - -^= (diV 3 (diipo, di 2 y 0 ) + d 2 V 3 {d 12 y 0 , d 22 y 0 )j ■ 
Using the relations between (dijyo , dj-yo) and diG in (|6.5I) . we obtain: 


(Vy 0 ) t Vp = 


1 

2A^2 


G 2 x 2 VU 3 < 8 > VA 


^=g 2 x 2 v 2 v 3 

V A 


1 

<vu 3 , VA) 

(VU 3 , VA X ) 1 

2\/A 

L — (VU 3 , VA 2 -) 

(W 3 ,VA) 
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and therefore: 

sym(V?/o) t Vp = y/X sym(VI/ 3 <8> V/) - \/AV 2 y 3 - V\{W 3 , VA)Id 2 . 

In a similar manner, it follows that: 

sym(Vy 0 ) 4 V4 = - A (v 2 / + |V/| 2 Id 2 ) . 

Since Q 2j a(z') = A _1 Q 2 , the formula in (13.511 becomes: 

WS) =h ^ e ' 2/ Q 2 ( s + ^(v^)*vy + f <g>v/) dz' 

+ ^7 [ e~ 2f Q 2 (2e / VR 3 <g> V/ - e f V 2 V 3 - e f {VV 3 , V/)Id 2 ) dx' 
24 Jq 

+ J n e_2/ Q 2 (e 2/ Ric(G) 2x2 ) dx', 


(7.3) 


which implies the result. ■ 

8. Appendix: a proof of Corollary 12.31 

1. For every x' € D denote D x .g = B(x',5 ) FI D and B x t^h = Ar',5 x (~h/2,h/2). For 
short, we write B^^h = B x ’. 2 h.h and B x ij t = B^ h ^. Apply Lemma [272] to the set V h = B x / }2 h 
to get a rotation R x ’ )2 h € SO( 3) such that, with a universal constant C: 


( 8 . 1 ) 


1 

h 


[ Vu h (z) — R x ', 2h (Qo(z') + z 3 Bo(z')) 


dz 


< C {^E h {u h ,B x i 2h) + h 3 |R x ', 2/ i|) • 

Consider a family of modifiers r\ x i G C°°( parametrized by x' G f2: 
r i (j (j 

J3 X ' = ^ II V x'%'IIl-(Q) < ^4 and (supp7/ x /) n n C 

Define € W 1)2 (D, M 3x3 ) as: 

(8.2) RV) = f r] x '(z')Vu h (z) (Qo(z') + zMz'))- 1 dz. 

Jn h 


We then have: 
1 
h 


[ | Vu h (z)-R h (z') (Qo(z') + z 3 B 0 (z')) I 2 dz 

JB. h 


< 


T f V/( 2 ) - R x ,, 2h (Q 0 (z') + z 3 B 0 (z')) 
h J B x ,, 2h 

+ % [ | R h (z') - R x ',2h\ 2 \Qo(z') + Z 3 B 0 (z ')\ 2 dz 

h JB. h 


,h 

< C (. E h (u h ,B x , !2h ) + h 3 |fl x ',2h|) +j l f B l#V) - 


(8.3) 
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where we have used (18.ip and \\Qo(z') + z 3 Bq(z')\\ Loo < C. Now, for every z' € B x > ^ we have: 


- R 


x',2h 


\R h {z')-R x ,, 2h \ 2 = [ r, z ,(y , )\7u h (y)(Q 0 (y , ) + y 3 B 0 (y')) 1 dy 
Jn h 

= J Vz'(y') (Vi/*(y) - R x ',2h (Qo(y') + y 3 -Bo(y'))) (QoOO + 2/3-60(2')) _1 d y 


( 8 - 4 ) < c 


riz'iy') dy 


IB 


z,h 



B 


Vu h {y)-R xl , 2h (Q 0 {y') + y 3 B Q (y')) 


z,h 


d y 


< S f | W*(y) - R x \ 2 h (Qo(y') + 2/360 (y')) 

J B x’, 2 h 

< Y 2 (E h (u h ,B xl>2h ) + h 3 \B x ,, 2h \) . 


d y 


In a similar way, in view of f Qh V z 'i] z r(y') d y = 0, it follows that: 

I VR h (z')\ 2 = U V z .r,Ai/)Vu h (y) ( Q 0 (y ') + yMy’))' 1 dy) 

V^y 2 /(y') (Vu h (y) (Qo(y') + y3-B 0 (y')) _1 - Rx', 2 h) dy 


I B 


< c f \V z >y z >(y')\ 2 dy [ \7u h \y) - R xJ , 2h (Q 0 (y') + y 3 B 0 (y')) 

B x',2h 

<^(E h (u h ,B x , t2h ) + h 3 \B x , }2h \y 

From (18.41) we obtain: 

[ \R h (z') - R x ,, 2h | 2 dz < ^2 [ (E\u\ B x , )2h ) + h 4 \B x ’ t2h \) dz 

J B x/ , h h JB x , h v y 

<C/l ( J F h (^,^, 2 / l ) + /i 3 |5x',2h|) , 

and therefore by (18.31) we further see that: 


dy 


(8.5) 


h j b 


x' ,h 


\Vu h (z)-R h (z') (Qo(z') + * 360 ( 2 ')) | 2 dz 

<c(E h (u h ,B x , !2h ) + h 3 \B xlj2h \) 


2. Covering fV' by a finite family of sets { B x > such that the intersection number of the 
doubled covering {B x / 2 h} is independent of h, applying (|8.5p and summing over the covering, 
it follows that: 

i J^ | Vu h (z) - R h (z') (Q 0 (zf) + z 3 6 0 (z')) I 2 dz < C ( E h (u h ) + h^ 
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In a similar fashion we obtain: 


VR h (z')\ 2 dz< 


Id 


x,h 


C_ 

hA 

C 


id. 


E h (u h , B x > t 2 h) + h 6 \B x ' 2 h \) dz 


< j~2 ( Eh (u h ,B x ,, 2h ) + h 3 \B x , j2h \) , 


and by the same covering argument: 


/ 


C 


Vi? ft (z')| 2 d z<A-( E h (u n ) + h 4 

nh h 2 


3. Note that, in the above two estimates, we can replace R h by R h = Pso( 3 ) R h € 
W 1,2 (Q., 50(3)). Firstly, the projection in question is well defined in view of (18.41) . since: 

dist 2 ( R h , 50(3)) < | R h - R x , t2h | < ^ (E h (u h ) + h 4 ) , 
which is small because of the hypothesis a < 2. Moreover: 


1 

h 



\Vu h (z) - R h {z ') (Qo(z') + z 3 B 0 (z’)) I 2 dz 

< t [ Vu h {z) - R h (z') (Qo(z') + z 3 B 0 (z')) 
h JB. h 


dz 


+ 


c 


L 


R h (z') - R h {z')\ 2 \Q G {z') + z 3 B 0 (z')\ 2 dz 


G|2 


< c(E h (u h ,B x , i2h ) + tf\B x ,, 2h \) 


because of ()8.5I) and (18.41) . Finally, the previous covering argument clearly implies (|2.9I) . and 
f n \VR h \ 2 dz < C \VR h \ 2 dz yields (107)1) . 
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